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We examine the properties of the Bethe Ansatz solvable two- and three-leg spin-S ladders. These 
models include Heisenberg rung interactions of arbitrary strength and thus capture the physics 
of the spin-S Heisenberg ladders for strong rung coupling. The discrete values derived for the 
magnetization plateaux are seen to fit with the general prediction based on the Lieb-Schultz-Mattis 
theorem. We examine the magnetic phase diagram of the spin-1 ladder in detail and find an 
extended magnetization plateau at the fractional value (Af ) = ^ in agreement with the experimental 
observation for the spin-1 ladder compound BIP-TENO. 
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I. INTRODUCTION 

After some intensive effort, the properties of the spin-i Heisenberg ladder are now well understood. The physics 
of the model is driven by the rung basis states, which dominate in the limit of strong rung coupling J . The spin gap 
originating between the singlet and triplet rung states for large J is known to persist for all finite J, vanishing only 
at the decoupling point J = where the model reduces to two Heisenberg chains. More generally, the n-leg spin-^ 
Heisenberg ladders possess the rather striking property of being massive for n even and massless for n odd. The 
ladder models also possess interesting magnetic properties, with clear evidence of magnetization plateaux. "^"^ From 
quite general arguments, the magnetization M of the spin-S* n-leg ladders should satisfy the discrete values 

nS{\ - M) e Z. (1) 

Although the Heisenberg ladders are not solvable in the sense of the spin-i Heisenberg chain, a number of solvable 
ladder models have been found. These include modified models with dimer or matrix product ground states^^^ and 
models solved by means of the Bethe Ansatz. Here we consider the latter class. These models include a two-leg 
spin-i model, ^ which compared to the pure Heisenberg ladder, has an additional four-spin interaction. The model 
is massless in the absence of rung interactions with a transition to a massive phase at a critical (non-zero) value of 
the rung coupling. In this solvable model the Heisenberg rung interactions appear as chemical potentials which break 
the underlying sm(4) symmetry. The phase diagram has been established by means of the Bethe Ansatz solution. ^'^^ 
In particular, magnetization plateaux are found at the values m = and m = 1, in agreement with (1). However, 
the Bethe Ansatz calculations also reveal an infinitcsimally thin magnetization plateaux at m = i, which may open 
up with the introduction of some anisotropy into the system. This thin plateau may be interpreted as a 'plateau 
boundary'. 

The solvable two-leg ladder model has been generalized to an arbitrary number n of legs in which the underlying 
su(2"') symmetry is broken by the Heisenberg rung interactions, which again appear as chemical potentials in the 
Bethe Ansatz solution. ^^■"'^^ Solvable ladder models with Hubbard"'^^ and t — J^^ interactions have also been found. 
Yet others have been constructed with mixed spin^^ and higher symmetries arising from other known integrable 
models. ^^'^^ Other models^* have been solved via mapping from the Hubbard model with appropriate boundary 
conditions. The su(4) ladder model has also been extended to include an additional free parameter. For the solvable 
3-leg spin-i model, ^^'^^ magnetization plateaux are found at the values m = ^ and m = 1 in agreement with (1).^*^ 
An additional infinitesimally thin plateau is also seen at m = ^. 

In this paper we investigate the properties of the Bethe Ansatz solvable two- and three-leg ladder models for general 
spin-S". These models are defined in Section II. In Sections III and IV we look respectively at the magnetization 
plateaux and their boundaries, and in Section V we give a detailed analysis of the two-leg spin-1 ladder. Our work on 
the solvable spin-1 ladder is motivated by the recent experimental and theoretical work on certain organic materials. 
For example, the compound BIP-TENO was studied experimentally,^^ and found to have magnetization plateaux of 
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0, ^, and 1. The crystal structure of this substance may be modeUed effectively by a two-leg spin-1 ladder. Some 
models have been used to propose mechanisms by which the fractional plateau appears. Relatively little work has 
been done on the spin-1 Heisenberg ladder. It has been studied in the semiclassical limit via a mapping to the 
nonlinear sigma model and by quantum Monte Carlo simulation and bosonisation techniques. ^^'^^ 
Concluding remarks are given in section VI. 



II. THE MODELS 



We consider the integrable n-leg spin-S ladders, or spin tubes with periodic boundary conditions. As usual a 'tube' 
is a ladder with an additional interaction between leg n and leg 1. We restrict our attention to the cases n = 2 and 
n — 3 with uniform spins on each leg, and not the more general mixed spin case. The spins along each leg and rung 
have an isotropic Heisenberg interaction, with the introduction of many-body terms to retain integrability. The su{2) 
spin-5 operators S = (5^, S^, S^) are defined by 

5^ = i(5+ + 5-), Sy = -li{S+-S-), 5^ = diag{5,5-l,...,-S}, (2) 

where, for given S, 



^m' m ~ V^iS + 1) - m{m - 1) S,n',m-l, (3) 



S-,,^ = VS{STl)^Mrr^+^) ^m',rn+l- (4) 

-S, —S +1,. . . ,S. The Hamiltonian of the model is 

H=f:Hij^i+j:Hr', (5) 

i=l i=l 
where L is the number of rungs. The leg term is given by 

H^L, = ® pS^ (n factors), (6) 

where 

n a-L 

a=0 PT^a 

with Xa = \oi{a + 1) — S{S + 1). Throughout this paper, site i + 1 is identified with site 1. For the simplest case of 
S = \, one recovers the familiar Heisenberg operator. Similarly for 5 = 1 and 5=1, 

Pf^, = {S,-S,+,f + S,-S,+,-l, (8) 
Pif+i = §('S'i • -Sj+i)^ -I- • - l^i • Sj+i - 11 . (9) 

The rung term has pure Heisenberg interactions, given by 



= E^'K-^r'-i)- (10) 

1=1 



where e = for tubes and e = 1 for ladders, and s\"'~^^^ = sf'K We take periodic boundary conditions along the 
ladder, and consider isotropic rung interactions {Ji = J). It has been shown that^^ + iJ™"^] = 0, so H 

defined in (5) is also integrable for sufficiently small n. 

A key ingredient in the models under consideration is that they are built from permutators. Such models can be 
diagonalized using the (nested) Bethe Ansatz.^^ The su{N) Bethe Ansatz equations are 

f^JljiV = |^ Ay)-Ai|)-i j| Af)-Ai^) + li ^ 
\A(^^ + Iv k=ij '^j^^ ~ + 1 fe=i '^f^ ~ ^^k^ ~ h 

Mr Ar) Jr) . M,_i ,(r) ,(r-l) i.M,+ i,(r) .{r+1) 1. 

n ^i ~ ^k -1 _ TT - ^k - 2^ TT ^3 - ^k - 2^ 

\W xM, • li xM x('-l) , 1: 11 xM , 1;' ^ ^ 

kjij ^/c "T 1 k=l ^fe "T 2^ k=l ^k ^ 2^ 



3 



where in this case N = {2S + 1)". As usual the eigenenergies of X^^i ^l^t+i 



Ml 



^ = ^-E77;^T2 — • (12) 



It is important to note that since the permutator (6) is invariant under any ordering of the basis states, this result 
may be obtained using any choice of reference state (or pscudo- vacuum) and any assignment of Bethe Ansatz 
pscudo particles. For each choice, however, one has to re-interpret the numbers Mr in terms of those corresponding 
to the ordering chosen. The rung Hamiltonian is dependent on the choice of ordering, but the change is just a 
rearrangement of its eigenvalues along the diagonal. We use this property to our advantage by doing calculations with 
that choice of ordering for which the Bethe Ansatz reference state is closest to the true groundstate of the system. 



A. Two-leg ladder 

We turn our attention now to the specific ease of the two-leg spin-S ladder, for the moment in the absence of a 
magnetic field. The rung Hamiltonian (10) is given by 

^rung ^ . g{2) _ ^ ^^3) 

which can be rewritten as 



^rung ^ J 



(g(l)+5(2))2_(g(l))2_(g(2))2 ^ 

2 



(14) 



It is convenient to change to the basis where the square and the ^-component of the total spin of a given rung are 

diagonal. Now, the operators {S^^^ + 5^^^)^, S^^'' , and S^"^^ are all diagonal with common eigenstates. As a result we 
can derive the eigenvalues of if ''"'^s simply by combining the eigenvalues of these operators as in the above expression. 
Prom elementary quantum mechanics, the respective eigenvalues of these operators are T(T -|- 1), S{S + 1), and 
S{S + 1), where the total spin T can assume the values 0, . . . , 25*. It follows that, for a two-leg ladder, the (25 -|- 1)^ 
rung states fall into 25-1-1 multiplets of total spin T, with eigenvalues 

Xt = ^{T{T+1)-2S{S+1)-2), T = 0,...,25, (15) 

and respective miiltiplicities rriT = 2T + 1 . 

It is emphasized that the Hamiltonian (6) has the same form in this new basis. The rung Hamiltonian (10) now 
becomes diagonal and is given by 

if™"s = diag{Ao,Ai,Ai,Ai,...,AT,...,AT,...,A2s,...,A2s}. (16) 

2T+1 4S+1 

In section V we will look at the specific example of a two-leg spin-1 ladder in detail. 

B. Three-leg tube 

We follow a similar procedure to the two-leg ladder, rewriting the Hamiltonian and changing the basis appropriately. 
In the case of the three-tube, the (25 -|- 1)^ states fall into [35J + 1 types of multiplet ([a;J denotes the integer part 
of x), with total spin T, whose eigenvalues are given by 

AT = ^(T(r+l)-35(5+l)-6), (17) 



where 



T='-^i^,'-^S-^ + l,...,3S. (18) 
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This has clear similarity to (15). However the situation now is somewhat more complicated. This is because in 
general there are several multiplets of given spin T; for example two distinct doublets having the same energy in the 
absence of a magnetic field. What follows is an analysis of exactly how the states fall into multiplets in a three-spin 
system. 

Let us consider the addition of three spin-S' objects. We first add the spins at sites (1) and (3) to get a combined 
spin-^ia system, and in turn, add the spin at site (2) to this combined system to got a state; of total spin T . Note 
that the choice of ordering of addition affects only the basis vectors used to represent a state, and not the actual 
state itself. The order chosen here facilitates computation of ladder eigenvalues in the next subsection. According to 
quantum mechanical rules for addition of angular momenta, we have 

< 5i3 < 25, (19) 

and 



|5i3-5| <r <5i3 + 5. (20) 

So the number of spin-T multiplets is given by the number of possible values of which satisfy (19) and (20). There 
are two cases to consider: 

case 1 : T < S 

In this case the right inequality of (20) is automatically satisfied. The left inequality is true when ^is assumes 
values in the range defined by (19) that satisfy 

|5i3 -S\<T, (21) 

namely 

Si3 = S-T,S-T+l,...,S + T; (22) 

and there are 2T + 1 such values. 
case 2: S <T 

We note that (19) implies that 1513 — S\ < S, so the left inequality of (20) is automatically satisfied. To make 
the right inequality true, 5i3 must assume values in the range defined by (19) that satisfy 

T < 5i3 + S, (23) 

namely 

5i3 = T-5, T-5+1, 25; (24) 
and there are 35 - T + 1 such values (since 25 = T - 5 + (35 - T)). 

To summarize these results: In a system of three spin-5 objects, the number of spin-T multiplets (that is, (2T -|- 1)- 
plets) is given by 

e(5,T) = | 5<T<35 (^^^ 

Multiplets of the same type have the same energy under the tube Hamiltonian, but are distinguishable by their value 
of 5i3, the possible values of which are 



'13 



5-r, 5-r + i,...,5 + r s>t , . 

T-S,T-S+1,...,2S S<T<3S. 



So, to compute the tube eigenvalues we use equation (17), and then use (25) to determine the number of multiplets 
which will have the given eigenvalues. 
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multiplet 


T 


quantity 




Sis from 


A«Si3 


shifts 


^{ladder) 


type 




(6(5, T)) 




(26) 








singlet 





1 


-6 


1 


-1 


2 


-4 


triplet 


1 


3 


-5 


0,1,2 


-2,-1,1 


3,2,0 


-2,-3,-5 


quintuplet 


2 


2 


-3 


1,2 


-1,1 


2,0 


-1,-3 


scptuplct 


3 


1 





2 


1 









TABLE I: Tabulated results relating tube eigenvalues and ladder eigenvalues. 



C. Three-leg ladder 

The thrcc-lcg ladder presents additional difRculty because of the broken symmetry. However, from (10) we see that 
the ladder and tube rung Hamiltonians are related via 

j^(ladder,rung) _ yj(tube,rung) _|_ j^^ _ ^ ^(3)-| ^27) 

Since S'^^^ ■ S^^^ commutes with the Hamiltonians (see Appendix A), their eigenvalues are similarly related via 

AS''^=Ar"^) + J(l-M5.3) (28) 

where usi^ is the eigenvalue of S^^^ ■ S^^^ for the state. Up to a shift and scaling, this operator is similar to the two-leg 
ladder Hamiltonian. By symmetry, the eigenvalues of S^^^ ■ S^^^ must be the same as those of S^^) • S'^'^\ and hence 
adaptation of (15) yields 

MSi3 = ^5i3(5i3 + 1) - S{S + 1) 5i3 = 0, 1, . . . , 2S (29) 

Thus, modifying the Hamiltonian from the tube to the ladder gives rise to an energy shift ( J(l — /isia)) ^^^^ depends, 
ultimately, on the value of S13 for the eigenstate (the tube and ladder Hamiltonians have the same eigenstates). So, 
to calculate the ladder eigenvalues, we proceed as follows: first calculate the tube eigenvalues along with the number 
of different type of multiplets that will occur using (25). Then, using (26) compute the range of values of 5*13 that 
will occur for a given type of multiplet. From these, use (29) to calculate fiSis ^^'^ hence the energy shift, and finally, 
add these to the original tube eigenvalues to obtain the ladder eigenvalues. 

The process is perhaps better clarified through an example, and we demonstrate for the case of the three leg spin-1 
ladder. We put J = 1 for simplicity. From (17) we determine that the tube eigenvalues are —6, —5, —3, and 0, with 
respective multiplicities 1 (one singlet), 9 (three triplets), 10 (two quintuplets), and 7 (a septuplet). The eigenvalues 
/txsi3 Ihis case are —2,-1, and 1, so the corresponding shifts (^(l — fJ^Sis)) required to give the tube eigenvalues are 
3, 2, and 0. 

We calculate the possible 6*13 values (and hence the possible ^Si3 values) that occur for each type of multiplet as 
given by equation (26), and then proceed to compute the energy shifts, and finally the ladder eigenvalues. The results 
are summarized in Table I. 



III. MAGNETIZATION PLATEAUX 



We now analyze the effect of introducing a magnetic field h along the positive ^-direction. This causes the multiplets 
to split into individual nondegenerate states according to the 2;-component of a state's spin. The relevant physical 

effect is the existence of magnetization plateaux, which occur when, in some region of phase space (defined by J and 
h), the groundstate consists of rungs all having the same magnetization. We work with the magnetization per site, 
which is defined as 

i=l i=l 

This definition is constructed in such a way that M has a saturation value of 1. As we move through phase space, 

passing through regions where the system is gapless, we find that the magnetization varies continuously. However in 
regions where there is a nonzero gap, the magnetization remains constant, hence the term 'magnetization plateau'. 
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The Lieb-Schultz-Mattis theorem for low- dimensional magnets gives rise to the necessary general condition^'^^ 

Q{Sc - rric) e Z, (31) 

for the existence of a plateau. Here Q is the spatial period as measured by the unit cell, and Sc and rric are the total 
spin and magnetization of the unit cell respectively. Hence for an n-leg ladder they may be taken as nS and nm. In 
this definition, m is not designed to have a saturation value of 1, and so we may take m = SM. For translationally 
invariant ladders, including those of this paper, we take Q = I, hence retrieving equation (1). 

In the case of the two-leg ladder, we note that the introduction of a magnetic field splits the energy levels of the 
states within a multiplet. Each state is shifted from the value given in (15) by an amount hS^, where is the 
z-component of its spin. Since the multiplet states may be characterized by their value of , there is a complete loss 
of degeneracy upon the introduction of a magnetic field. For a rung of total spin T, the lowest energy state from each 
one has eigenvalue 

^T,imn = ^iT{T+l)-2S{S+l)-6]-hT, T = 0, . . . , 25. (32) 

We show that it is possible for any of these states to be the groundstate with a nonzero gap, subject to appropriate 
choice of J and h. Consider first the case where the two rung states of lowest energy are from the same multiplet. 
Since the magnetic energy shifts are given by /i5^, it is clear that the difference in their energy due to the rung 
Hamiltonian is h. The leg Hamiltonian gives a contribution of 4, so that the energy gap is A = h — 4. Therefore if 
h> 4, all rungs are in the same state, and the system has a gap A. 

In the case where the two lowest energy states are from diff'erent multiplets (say, Ap and Xq, where Aj, > Xq), the 
gap is given by 

A = Ap-A, -4 

= |b(p+i)-a(« + i)]-%-a)-4. (33) 

This can be guaranteed to be greater than zero for any p and q simply by taking 

[p{p+l)-q{q+l) 

in the antiferromagnetic case (with a similar result in the ferromagnetic case). Hence, the lowest energy state from each 
multiplet has some corresponding region of phase space where it is the sole rung state occurring in the groundstate, 
and the ladder is gapped. The different multiplets take on all possible values of Sz allowed under the restriction that 
the system consists of intrinsic spins 5. Consequently, it is apparent that a plateau of every allowable magnetization 
is realized somewhere in phase space. This analysis can also be applied to the three-leg case without significant 
modification. 

It is simple to check that the magnetization plateax we have found here satisfy the requirements of (31). We note 
that, if the gaps discussed above are nonzero, then all rungs are in the same groundstate. Hence Q = 1, and we 
require that S — m ^ k, where 5 and m are the spin and magnetization of a rung, and k is an integer. However, since 
m can take values ranging from —5 to S in integer steps, it is clear that this relation will be satisfied. 



(34) 



IV. PLATEAU BOUNDARIES 



It is of some interest to explore what happens on the boundary between two magnetization plateaux. We examine 
the two-leg case. The boundary may be defined by the line on which the groundstate of each plateau is degenerate. 
In general, we see that the states that can have minimum energy in the quadrant J,h > are the states from each 
multiplet with largest z-component of the spin, namely 



^T,min = ^[TiT+l) - 25(5+1)] - hT. (35) 

It is also clear that as we steadily increase the ratio h/J that T will progress in order through its possible values 
0, .., 25. That is, the plateau boundaries will occur between successively ordered states as in (35). Thus the equations 
of the plateau boundaries are given by Ar.min — Ar+i.min = or more explicitly 



h=-{T+l) T = 0,l,...,25-1. 



(36) 
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Now it is a simple matter to determine what the magnetization wih be on these hnes. As usual we normalize so 
that {M) = 1 if 5^ = 5* on each site. We calculate the magnetization on these boundary lines by averaging the two 
states that meet there. Now the state corresponding to AT.min must have Sz as large as possible, that is, Sz = T. On 
boundary line T, then, we simply average T and T + 1, which is spread over two sites and normalized for spin. That 
is 

\ )/2S = T = 0,l,...,25-1. (37) 

Although we would not, strictly speaking, consider this a magnetization plateau, it still satisfies (31) in some sense. 
We can think of the groundstate now as being a mixture of two types of rung. Hence for the two-leg ladder, we take 
Q = 2 to be a unit cell. The magnetization of the unit cell is then 2T + 1, with the total spin being some integer 
between and 45. Hence 2{S — m) as defined in (31) is clearly an integer. 



V. EXAMPLE - THE TWO-LEG, SPIN-1 LADDER 



We look at the specific case of the two-leg, spin-1 ladder in the presence of a magnetic field. In this example the 
leg Hamiltonian is 



5W.5W + (5W.5«)'-1 



^(2) ^ ^(2) 



while the rung Hamiltonian is given by 

= J (Sf^ . -l)-h ({S^f^ + {S^)f 

where 

/O 1 o\ 



(38) 



(39) 



1 

71 



1 1 

Vo 1 oy 



sy = 



72 



/O -i 0\ 




(I 





o\ 


i -i 













Vo i oj 




lo 





-1/ 



(40) 



are the spin-1 operators. We will work in the basis where the square and the ^;-component of the total spin of a given 
rung, S = 5^^^ -I- 5^^^ , are diagonal. It follows that the nine states on a given rung fall into a spin-2 quintuplet, a 
spin-1 triplet, and a spin-0 singlet. Let S = \S\ and let be the z-component of the total spin. We denote the rung 
basis by \S, S^) and it is given in terms of the (5^)^'^ eigenstates by 



|0) = 


|0; 


0) = 


75(l+)|-)-|0)|0)-f 1- 




ll) = 


|1; 


1> = 


;^(-|+)|0) + |0)|+)) 




|2) = 


|1; 


o> = 






|3) = 


|1; 


-1> 


= ;^(-|0)|-) + |-)|0)) 




|4) = 


|2; 


2) = 


l+>l+) 




|5) = 


|2; 


1) = 


^(l+)|0) + |0)|+)) 




|6) = 


|2; 


0) = 


^(|+)|-)+2|0)|0) + | 




|7) = 


|2; 


-1) 


= ;^(|0)|-) + |-)|0)) 




|8) = 


|2; 


-2) 


= |-)|-) 





(41) 

The leg Hamiltonian (6) has the same form in this new basis. The rung Hamiltonian (39) now becomes diagonal and 
is given by iJ™"*? = diag{-3 J, -2 J - h, -2 J, -2 J + h, -2h, -h, 0, h, 2h}. 

We now consider the possible orderings of the states according to their energy in the presence of the magnetic field h. 
The energy considerations in this case are quite simple, as there is no degeneracy. The number of possible orderings is 
now quite large because the field can cause varying amounts of energy level splitting within a given multiplet, leading 
to overlaps between them. In fact there are 12 physically possible orderings for h > Q. From now on we will restrict 
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our attention to the first quadrant (J > 0) which contains six orderings. We will list these below, computing energy 
gaps, and hence deriving the phase diagram. 

In computing energy gaps for each ordering, we will use a procedure similar to that used in work on the su{8) 
tube.^^ Since the rung Hamiltonian is diagonal, it is a simple matter to compute its contribution to this gap. However, 
calculating the contribution of the leg Hamiltonian is more complicated. This requires taking the thermodynamic 
limit of the Bethe Ansatz ciquations, taking Fourier transforms, and then solving. The final step comes down to the 
evaluation of an integral, which gave us our final result. 

This integral, and hence the energy gap, ultimately depends on the type of mode we need to create in order to 
compute the energy gap. This follows directly from the degeneracy of the groundstate. It is therefore possible to give 
a general form of this integral, which depends only on the type of mode we are creating. This allows us to easily and 
quickly compute the gap contribution of a Hamiltonian that is to be solved by the Bethe Ansatz method. 

If we denote by afe the integral corresponding to creating a A^*^^ excitation, the general result is 



-2 

Oik = -r 

k 



7 + log4 + V(i + ^) 



(42) 



where 7 is the Euler gamma constant and ^(2^) is the digamma function. This can be simplified for a specified k. 
Here, we will only need the values ai = 4 and a2 = 2 log 2 = 1.38629 . . . 
We now proceed to analyze each ordering in turn. 

1. h<lJ 

In this case the ordering is {|0), |2), |3), |4), |5), |6), |7), |8)}. The energy is given by = Ei^g + Mi{J - h) + 
(M2 + Ms)/! + Ma{2J - 3h) + (M5 + Me + M7 + Ms)h, where Mk = J2l=k ^Pr^ where A'' is the number of 
states |j) ocurring in an excited state, and P takes the numbers 0, 1, 2, 3, 4, 5, 6, 7, 8 to the ordering we present. 
(In this case, P is the identity permutation, in the cases below, it will be some other permutation). Hence 
the first energy gap is A = (— 2J — h) — (— 3J) — 4 = J — h — A, leading to a critical line he = J — 4:. There 
is another crosssover curve on which the second excitations, (|2) states) also become massless. Although we 
cannot calculate this exactly, we can see that it intersects the above critical line at the point {J,h) = (4,0). 
This curve moves as Mi changes, eventually traversing regions 1. and 2. to the point {J,h) = (2 log 2, 2 log 2) 
on the boundary of region 3, h = J. 

Below the line /ic = J — 4, all rungs are in state |0), so the total magnetization of the system is (M) = 0, as 
would be expected. Above this line a finite number of rungs exist in the lowest-energy triplet state, and the 
magnetization varies continuously. 

2. IJ <h<J 

Now the ordering is {|0), |1), |2), |4), |3), |5), |6), |7), |8)}. This gives the same calculation as for case 1, so the 
critical line is the same. Regions 1. and 2. become distinguishable only when considering the third or higher 
excitations. 

3. J < /i< |J 

Now the ordering is |0), |4), |2), |5), |3), |6), |7), |8)}. The first energy gap is .B = - J + /i - 4, giving rise to 

a critical line of /ic = J + 4. The second excitations become massless on this line at (J, h) = (8, 12), and move 
to the point {J,h) = (2 log 2, 2 log 2). That is, the same point as that calculated for region 1. 

Now above this line, the groundstate consists entirely of rungs in state |1). As a result there is a magnetization 
plateau of (M) = ^ here. 

4. p <h<2J 

The ordering is {|1), |4), |0), |2), |5), |3), |6), |7), |8)}. The gap now is A = 2J - ft - 4 and the critical curve is 
/ic = 2 J — 4. This intersects the critical curve of region 3. on the boundary between them. Again the second 
excitations become massless at (J, /i) = (8,12), and move to the point {J,h) = (2 log 2, 4 log 2) on the line 
h = 2J. 

Below the line he = 2 J — 4 the groundstate consists entirely of rungs in the |1) state, yielding a magnetization 
of (M) = i, being part of the same plateau as found in region 3. 

5. 2J < h <3J 

We now have an ordering of {|4), |1), |0), |5), |2), |6), |3), |7), |8)}. The critical line here is he = 2J + 4. The 
second excitations do not become massless on this curve. They do, however, become massless at the point 
(J,/i) = (21og2,41og2). 
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FIG. 1: The magnetic phase diagram of the two-leg spin-1 ladder. The thin lines are h = ^J, h = J, h = ^J, h = 2J and 
h — 3J, which together with the y-axis divide the quadrant into the regions 1-6. The thick lines are the exactly derived phase 
boundaries, and the curved lines are sketches of higher phase boundaries. 

6. /i > 3 J 

The ordering is {|4), |1), |5), |0), |2), |6), |3), |7), |8)}, yielding an extension of the critical line found in case 5. The 
second excitations become massless on this line at (J, h) = (0, 4). 

Above the line = 2J + 4 in regions 4. and 5. the groundstate consists entirely of rungs in state |4), and the 
system reaches a magnetization plateau of (M) = 1. Outside the three plateaux found here the magnetization 
varies continuously. 

The results of this analysis are encapsulated in the phase diagram (Fig. 1). We show its first quadrant. The 

shape of the curved line is interesting, showing a sharp spike with its cusp on the line h = | J. Although we have 
calculated some exact points, its overall shape cannot be determined analytically. We can predict some features using 
perturbation theory, and it is possible that the shape of the curve could be accurately determined using numerical 
scries methods. Note also that there is more structure inside this curve, corresponding to the points at which higher 

excitations become massless. 



VI. CONCLUDING REMARKS 

We have investigated the behaviour of some Bethe Ansatz solvable two- and three-leg ladder models with general 
spin S. These models include multi-spin interactions in order to be exactly solvable. Most importantly, the solvable 
ladder models include Heisenberg rung interactions of arbitrary strength which drive the physics for strong rung 
coupling. The values of the magnetization on each of the magnetization plateaux are seen to fit with the general 
prediction based on application of the Lieb-Schultz-Mattis theorem. Additionally, the plateau boundaries, which may 
be thought of as infinitesimally thin plateaux, also satisfy this criterion. We calculated the phase diagram (Fig. 1) of 
the spin-1 ladder in detail. Distinct plateaux occur at the values (M) = 0, i, 1 in agreement with the experimentally 
observed^° values for the spin-1 ladder compound BIP-TENO. It will be most interesting to derive the thermodynamic 
properties of this model. 

This work has been supported by the Australian Research Council through grants A6990558, F69700124 and 
DP0208925. M. Maslen has been supported by an Australian Postgraduate Research Award. 
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APPENDIX A: COMMUTATION OF S^^' S^^^ WITH THE HAMILTONIANS 

We first show that S^^^ ■ S^^^ commutes with the rung Hamiltonian. For the three-leg case we have 

jjrung ^ _j(-2 _ . S(2) _ 5(2) . 5f(3)) ^ 



(Al) 



so 



= J 
= J 



5(1).5(3) ,5(l).;5(2)+5(2).^(3)j^ . 

This commutator may be expressed as a linear combination of commutators of the form 



(A2) 



(A3) 



where a, b and c represent the co-ordinate variables x, y and z, and k, I, m represent the site number. Using standard 
properties of commutators, this can be expanded to give 



■(^a)(fc) ^ ^^b^ Wl ^gaf) ^gb-ji"^) ^ If^gaf) ^ (^.c^lOl (^a^Cfe) ^gc^{m) 

and, using the quantum mechanical relation [Sj,Sk] = iejkiSi, this becomes 

i£a6c (5^)^'^ (S*-)^"^ + i£„e6 (5*') (5^)^™) 

= 0, 



(A4) 



(A5) 



where summation is implied over repeated indices, and we have used the fact that spins from different sites commute. 

That S*-^-* -S^"^^ commutes with the leg Hamiltonian follows directly from equation (17) of a previous paper. -"^^ As was 
noted after that equation, arbitrary XYZ interactions introduced on the rungs will commute with leg Hamiltonians 
of the type we are working with here. Of course, the same statement is true for removing such interactions, which 
is how the tube Hamiltonian is transformed into the ladder Hamiltonian. Since the Hamiltonians and the operator 
5(1) . 5(3) are mutually commutative, the eigenvalues may be related directly as in (28). 
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